Coherent dynamics of a flux qubit coupled to a harmonic oscillator 
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In the emerging field of quantum computation 1 
and quantum information, superconducting devices are 
promising candidates for the implementation of solid- 
state quantum bits or qubits. Single-qubit operations 2-6 , 
direct coupling between two qubits 7 ~ 10 , and the real- 
ization of a quantum gate 11 have been reported. How- 
ever, complex manipulation of entangled states — such 
as the coupling of a two-level system to a quantum 
harmonic oscillator, as demonstrated in ion/atom-trap 
experiments 12,13 or cavity quantum electrodynamics 14 — 
has yet to be achieved for superconducting devices. Here 
we demonstrate entanglement between a superconduct- 
ing flux qubit (a two-level system) and a superconducting 
quantum interference device (SQUID). The latter pro- 
vides the measurement system for detecting the quantum 
states; it is also an effective inductance that, in parallel 
with an external shunt capacitance, acts as a harmonic 
oscillator. We achieve generation and control of the en- 
tangled state by performing microwave spectroscopy and 
detecting the resultant Rabi oscillations of the coupled 
system. 

The device was realized by electron-beam lithogra- 
phy and metal evaporation. The qubit-SQUID geom- 
etry is shown in Fig. la: a large loop interrupted by 
two Josephson junctions (the SQUID) is merged with 
the smaller loop on the right-hand side comprising three 
in-line Josephson junctions (the flux qubit) 15 . By apply- 
ing a perpendicular external magnetic field, the qubit 
is biased around $o/2, where $o = h/2e is the flux 
quantum. Previous spectroscopy 16 and coherent time- 
domain experiments 6 have shown that the flux qubit is 
a controllable two- level system with 'spin- up/spin-down' 
states corresponding to persistent currents flowing in 
'clockwise/anticlockwise' directions and coupled by tun- 
neling. Here we show that a stronger qubit— SQUID cou- 
pling allows us to investigate the coupled dynamics of a 
'qubit— harmonic oscillator' system. 

The qubit Hamiltonian is defined by the charging and 
Josephson energy of the qubit outer junctions (Ec = 
e 2 /2C and Ej = hlc/ie where C and Ic are their capac- 
itance and critical current) 16 . In a two- level truncation, 
the Hamiltonian becomes H q /h = —e<j z /2—Ao- x /2 where 
a z%x are the Pauli matrices in the spin- up/spin-down ba- 
sis, A is the tunnel splitting and e = 7 p $o(7 9 — ir)/hn 
(Ip is the qubit maximum persistent current and 7 g is 
the superconductor phase across the three junctions). 



The resulting energy level spacing represents the qubit 
Larmor frequency Fl — \/A 2 + e 2 . The SQUID dynam- 
ics is characterized by the Josephson inductance of the 
junctions Lj ss 80 pH, shunt capacitance C s h ~ 12 pF 
(see Fig. la) and self-inductances L s i w 170 pH of the 
SQUID and shunt-lines. In our experiments, the SQUID 
circuit behaves like a harmonic oscillator described by 
H sq = hv p (ala + 1/2), where 2irv p = l/y/(Lj + L sl )C sh 
is called the plasma frequency and a (a T ) is the plasmon 
annihilation (creation) operator. Henceforth \(3n) rep- 
resents the state with the qubit in the ground (/3 = 0) 
or excited ((3 = 1) level, and the oscillator on the n th 
level (n = 0, 1, 2, . . .). The corresponding level diagram 
is sketched in Fig. lb (inset). The coupling between the 
qubit and the oscillator originates from the current dis- 
tribution in the shared branches (Fig. la) and gives rise 
to an interaction Hamiltonian H q - sq = Xu z (a + a') with 
A w 0.2 GHz in our device 17 (the estimated qubit-SQUID 
coupling is M w 20 pH). 

Measurements are performed at T=25 mK using low- 
noise circuitry to minimize decoherence, relaxation and 
thermal activation. The system is first initialized by al- 
lowing it to relax to the |00) ground state. With suc- 
cessive resonant microwave pulses we achieve controlled 
superposition of various \/3n) states, as shown below. 
The readout 6 is performed by applying a short current 
pulse lb (~30 ns) and by monitoring whether the SQUID 
switches to the finite-voltage state. After averaging ty- 
pically 10000 readouts, we obtain the probability P sw (Ib) 
which for properly-chosen parameters is proportional to 
the excited state occupancy. In the following, we first 
show the spectroscopy of the coupled qubit— oscillator 
system and Rabi oscillations of the qubit. Next we 
demonstrate coherent dynamics of the coupled system. 

We performed spectroscopy of the coupled 
qubit— oscillator system by applying a long (300 ns) 
microwave pulse with various frequencies and measuring 
the SQUID switching probability. Peaks and dips are 
observed and their resonant frequencies as a function of 
A$ = $ ext — $0/2 are given in Fig. lb. We obtain one 
manifold of three resonances spaced by ~2.91 GHz. This 
frequency coincides with the designed oscillator eigenfre- 
quency v v . In addition, we observe a spectroscopic peak 
or dip that depends only weakly on the magnetic field 
(circles in Fig. lb). For lower microwave power, only 
the qubit band (squares) remains visible. A numerical 
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FIG. 1: Qubit— SQUID device and spectroscopy a, Atomic 
force micrograph of the SQUID (large loop) merged with 
the flux qubit (the smallest loop closed by three junctions); 
the qubit to SQUID area ratio is 0.37. Scale bar, lfxm. 
The SQUID (qubit) junctions have a critical current of 4.2 
(0.45) fiA. The device is made of aluminium by two sym- 
metrically angled evaporations with an oxidation step in be- 
tween. The surrounding circuit shows aluminium shunt ca- 
pacitors and lines (in black) and gold quasiparticle traps 3 and 
resistive leads (in grey). The microwave field is provided by 
the shortcut of a coplanar waveguide (MW line) and couples 
inductively to the qubit. The current line (/) delivers the 
readout pulses, and the switching event is detected on the 
voltage line (V). b, Resonant frequencies indicated by peaks 
in the SQUID switching probability when a long microwave 
pulse excites the system before the readout pulse. Data are 
represented as a function of the external flux through the 
qubit area away from the qubit symmetry point. Inset, en- 
ergy levels of the qubit— oscillator system for some given bias 
point. The blue and red sidebands are shown by down- and 
up-triangles, respectively; continuous lines are obtained by 
adding 2.96 GHz and —2.90 GHz, respectively, to the central 
continuous line (numerical fit) . These values are close to the 
oscillator resonance u p at 2.91 GHz (solid circles) and we at- 
tribute the small differences to the slight dependence of v v on 
qubit state, c, The plasma resonance (circles) and the dis- 
tances between the qubit peak (here Fl = 6.4 GHz) and the 
red/blue (up/down triangles) sidebands as a function of an 
offset current hoff through the SQUID. The data are close 
to each other and agree well with the theoretical prediction 
for i/p versus offset current (dashed line). 




MW pulse length (ns) F, (GHz) 

FIG. 2: Rabi oscillations at the qubit symmetry point A = 
5.9 GHz. a, Switching probability as a function of the mi- 
crowave pulse length for three microwave nominal powers; de- 
cay times are of the order of 25 ns. For A = 8 dBm, bi-modal 
beatings are visible (the corresponding frequencies are shown 
by the filled squares in b). b, Rabi frequency, obtained by 
fast Fourier transformation of the corresponding oscillations, 
versus microwave amplitude. In the weak driving regime, the 
linear dependence is in agreement with estimations based on 
sample design. A first splitting appears when the Rabi fre- 
quency is ~v p . In the strong driving regime, the power in- 
dependent Larmor precession at frequency A gives rise to a 
second splitting, c, This last aspect is obtained in numer- 
ical simulations where the microwave driving is represented 
by a term [\/2)hF\ cos(Af) and a small deviation from the 
symmetry point (100 MHz) is introduced in the strong driv- 
ing regime (the thick line indicates the main Fourier peaks). 
Radiative shifts 20 at high microwave power could account for 
such a shift in the experiment. 

fit (continuous line) of this band leads to Ej = 225 GHz, 
Ec = 7.3 GHz, and ratio of area of qubit junctions 
a = 0.76 (A = 5.9 GHz, I p = 275 nA). The appearance 
of the manifold instead of a single resonance is due to 
the qubit coupling with the oscillator mode f p (ref. 18). 
Similarly to atomic physics, we call the |00) to |11) 
( 1 01) to 1 10)) transitions the blue (red) sidebands (see 
the ladder energy diagram of the \(3n) states in Fig. lb 
inset). We note that near the qubit symmetry point, 
the closeness of the oscillator resonance and the red 
sideband, visible owing to a small thermal occupation of 
the 1 01) state, is purely accidental. To verify that the 
oscillator involved is indeed the SQUID plasma mode, 
we repeated the above measurements in the presence of 
an offset bias current hoff which decreases the plasma 
frequency following 19 v p {\ — (Iboff/Ic) 2 ) 1 ^ 4 , where I c 
is the SQUID critical current (4.2 A). The data in 
Fig. lc show the distance between the qubit peak for 
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FIG. 3: Generation and control of entangled states, a, Spec- 
troscopic characterization of the energy levels (see Fig. lb 
inset) after a n (upper scan) and a 2tt (lower scan) Rabi pulse 
on the qubit transition. In the upper scan, the system is first 
excited to 1 10) from which it decays towards the 1 01} excited 
state (red sideband at 3.58 GHz) or towards the |00) ground 
state (Fl = 6.48 GHz). In the lower scan, the system is ro- 
tated back to the initial state 1 00) wherefrom it is excited into 
the 1 10) or 1 1 1} states (see, in dashed, the blue sideband peak 
at 9.48 GHz for 13 dB more power), b, Coupled Rabi oscilla- 
tions: the blue sideband is excited and the switching proba- 
bility is recorded as a function of the pulse length for different 
microwave powers (plots are shifted vertically for clarity). For 
large microwave powers, the resonance peak of the blue side- 
band is shifted to 9.15 GHz. When detuning the microwave 
excitation away from resonance, the Rabi oscillations become 
faster (bottom four curves) . These oscillations are suppressed 
by preparing the system in the 1 10) state with a 7r pulse and 
revived after a 2-k pulse (top two curves in Fig. 3b) c, Cou- 
pled Rabi oscillations: after a it pulse on the qubit resonance 
(| 00) 1 10)) we excite the red sideband at 3.58 GHz. The 
switching probability shows coherent oscillations between the 
states 1 10) and 1 01} , at various microwave powers (the curves 
are shifted vertically for clarity). The decay time of the co- 
herent oscillations in a, b is ~3 ns. 



Fl = 6.4 GHz and the blue/red sidebands (down/up 
triangles) that decreases together with the oscillator 
resonance (circles). 

To realize quantum operations on the qubit only, we 
apply a resonant microwave pulse with frequency F mw — 
Fl. The operation is performed at the qubit symmetry 
point -f q = n where F L = A. In Fig. 2a, the SQUID 
switching probability is plotted against the microwave 
pulse length for three microwave power levels. The ob- 




2nd pulse length (ns) delay 2-1 (ns) 



FIG. 4: Oscillator relaxation time, a, Rabi oscillations be- 
tween the |01) and 1 10) states (during pulse 2 in the inset) 
obtained after applying a first pulse (1) in resonance with 
the oscillator transition. Here, the interval between the two 
pulses is 1 ns. The continuous line represents a fit using an 
exponentially decaying sinusoidal oscillation plus an exponen- 
tial decay of the background (due to the relaxation into the 
ground state). The oscillation's decay time is r co h = 2.9 ns, 
whereas the background decay time is ~ 4 ns. b, The ampli- 
tude of Rabi oscillations as a function of the interval between 
the two pulses (the vertical bars represent standard error bars 
estimated from the fitting procedure, see a). Owing to the os- 
cillator relaxation, the amplitude decays in r re ; r; 6 ns (the 
continuous line represents an exponential fit). 

served Rabi oscillations decay within ^30 ns. Remark- 
ably, we can reach Rabi frequencies comparable to the 
Larmor frequency (up to 6.6 GHz). Using Fourier trans- 
formation, we extract the Rabi frequency as a function of 
the microwave amplitude (Fig. 2b). In the weak driving 
regime, the Rabi frequency increases linearly with the 
microwave power, as expected 6 . Near the oscillator res- 
onance v v , we see two frequencies in the spectrum, a be- 
haviour which is probably caused by the qubit— oscillator 
coupling. At even higher microwave powers, the spec- 
trum exhibits again a second frequency component at 
A. A qualitatively similar behaviour is also obtained 
in numerical simulations (see Fig. 2c) when we consider 
the qubit driven by an additional term (l/2)/iFi cos(Ai) 
in H q (Fi and A are the microwave amplitude and fre- 
quency, respectively). 

We now turn to the conditional dynamics resulting 
from the qubit— oscillator coupling. We first determine 
the blue and red sideband resonant frequencies by spec- 
troscopic means using a two-pulse sequence (Fig. 3a). 
The qubit is prepared in the excited state by a n pulse at 
the Larmor frequency. A second pulse (18 ns) of variable 
frequency induces resonant qubit de-excitation (dips in 
Fig. 3a top trace) marking the red sideband and the Lar- 
mor frequency. Similarly, after a 2-7T pulse which places 
the qubit in its ground state, we search for resonant ex- 
citations (peaks in Fig. 3a bottom trace) that mark the 
Larmor frequency and the blue sideband. No resonance 
is seen on the red sideband, showing that the oscillator 
is in its ground state with a large probability. Note that 
in order to excite the blue sideband, we have to increase 
the microwave power by at least 10 dB, probably due to 
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less effective microwave transmission in the 8 — 9 GHz 
range (note also the absence of spectroscopy peaks in 
this frequency range in Fig. lb). At high microwave 
powers, we observe radiative shifts 20 of the resonances. 
We now exploit these resonances to study the dynam- 
ics of the coupled system by applying pulses of varying 
length. In Fig. 3b, Rabi oscillations are shown for the 
|00) to 1 11) transition. When the microwave frequency 
is detuned from resonance, the Rabi oscillations are ac- 
celerated (bottom four curves, to be compared with the 
fifth curve). After a tt pulse which prepares the system 
in the 1 10) state, these oscillations are suppressed (sec- 
ond curve in Fig. 3b). After a 2tt pulse they are revived 
(first curve in Fig. 3b). In the case of Fig. 3c, the qubit 
is first excited onto the 1 10) state by a tt pulse and a sec- 
ond pulse in resonance with the red sideband transition 
drives the system between the 1 10} and 1 01} states. The 
Rabi frequency depends linearly on the microwave am- 
plitude, with a smaller slope compared to the bare qubit 
driving. During the time evolution of the coupled Rabi 
oscillations shown in Figs. 3b and 3c, the qubit and the 
oscillator experience a time-dependent entanglement, al- 
though the present data do not permit us to quantify it 
to a sufficient degree of confidence. 

The sideband Rabi oscillations of Fig. 3 show a short 
coherence time (^3 ns) which we attribute mostly to the 
oscillator relaxation. To determine its relaxation time, 
we performed the following experiment. First, we excite 
the oscillator with a resonant low power microwave pulse. 
After a variable delay At, during which the oscillator re- 
laxes towards n = 0, we start recording Rabi oscillations 
on the red sideband transition (see Fig. 4a for At = 1 ns). 
The decay of the oscillation amplitude as a function of 
At corresponds to an oscillator relaxation time of ^6 ns 
(Fig. 4b), consistent with a quality factor of 100 — 150 
estimated from the width of the v p resonance. The ex- 
ponential fit (continuous line in Fig. 4b) shows an offset 
of ^4% due to thermal effects. To estimate the higher 
bound of the sample temperature, we consider that the 
visibility of the oscillations presented here (Figs. 2-4) is 
set by the detection efficiency and not by the state prepa- 
ration. When related to the maximum signal of the qubit 
Rabi oscillations of ~40%, the 4%-offset corresponds to 
~10% thermal occupation of oscillator excited states (an 
effective temperature of ~60 mK). Consistently, we also 
observe low-amplitude red sideband oscillations without 
preliminary microwave excitation of the oscillator. 

We have demonstrated coherent dynamics of a cou- 
pled superconducting two-level plus harmonic oscillator 
system, implying that the two subsystems are entangled. 
Increasing the coupling strength and the oscillator relax- 
ation time should allow us to quantify the entanglement, 
as well as to study non-classical states of the oscillator. 
Our results provide strong indications that solid-state 
quantum devices could in future be used as elements for 
the manipulation of quantum information. 
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